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Abstract 

We study the decay properties of a heavy-light meson. We reformulate the decay 
amplitudes for the heavy-light systems and find a new way to calculate decay rates. 
Applying this formulation, we find a new sum rule for the radiative decays of one 
heavy-light meson into another, H\ — > H2 + 7 with various combinations of H^. 
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1. Introduction 



Ever since the discovery of the so-called D s j in 2003 by BaBar 11 and CLEO, 2 ' which are 
now identified as £> s (2317) and £> s (2460) with = (1/2)+ and (3/2)+, respectively, the 
particles are discovered in the vicinity of this energy or above one after another. Theorists 
as well as experimentalists are trying to explain these particles as ordinary Qq, or molecular, 
or tetra-quark, or a bound state of diquark states, eteP At the present time there are several 
explanations for the same state and the final decisive answer is not yet given. 

As well as explanations of its origin, many people study the properties of these heavy 
mesons, spin, parity, decay rates, etc. In this paper, we report our study on the formulation 
to calculate decay rates of a heavy-light meson which are Lorentz-invariant and on the kine- 
matical results derived from it. We assume that the D s j are heavy-light composite mesons 
Qq, whose assumption is supported by experimentally compatible numerical calculations of 
the lattice gauge theory.®'® 

There are several ways to calculate decay rates of the heavy-light mesons proposed so 
far; one is by Goity and Roberts® which is applied to newly discovered heavy mesons by Di 
Pierro and EichtenP and another by Bardeen, Eichten, and Hill. 8 ' The former is based on 
the ordinary dipole expansion of a bound state with the heavy quark at rest at the center 
and the latter utilizes the effective chiral Lagrangian which describes coupling of heavy-light 
mesons with chiral particles. 

We propose a different method from these to calculate decay rates, which are Lorentz- 
invariant and are calculated in the Breit frame where initial and final mesons have the same 
velocity in the opposite directions. 



We start to construct formulation how to calculate a decay rate of a bound state because 



obtain the numerical results. We would like to solve these problems in this paper. For 
instance, consider the method proposed in Refs. EJ and [7j), in which, e.g., pion decay is 
described by the following equation. 



§2. Formulation 



the former methods are not Lorentz-covariant®^ or do not use meson wave functions to 





where the interaction is assumed to be 



int 
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with axial vector current js M . Here one uses the composite meson wave functions, ipij, at the 
rest frame even for the final state ipf. Another point is that a plane wave pion wave function 
exp (ik ■ xj is inserted ad hoc. Without phase factors in the initial and final wave functions 
exp (—iPij ■ x), the expression of Eq. (12-11) neglects the recoil of heavy-light mesons. We 
claim the factor in Eq. (12 -lp should be replaced with the following expression in the Breit 
frame. 



exp (—ik ■ x) — > exp (— HrriQV z) 



(2-3) 



where tuq is a heavy quark mass. The factor 2 on the r.h.s. of Eq. (12-31) in the exponent 
appears because the initial heavy quark has a momentum rriQ^/V in the +z direction and 
the final one has rriQ'jV in the — z in the Breit frame, hence the recoil momentum becomes 
2rriQ r )V . The factor 7 is absorbed into integral variable z' = 72. 

Both sides of Eq. (I2-3P are equal to each other when one replaces uiq on the r.h.s. with an 
average of hadron masses as (mi + 777.2) /2. When we include phase factors of initial and final 
wave functions ipij, together with the plane wave pion wave function, we obtain four-delta 
function meaning four-momentum conservation. That is, the plane wave pion wave function 
disappears from the final expression but there remains the recoil phase factor. 

Our derivation of Eq. (I2-3|) is based on the field theory and is given as follows. 



(0\q c (t,x)Q(t,y)\P) 







q c (0,x-Xt)Q(Q,y-Xt 



P)e 



i.P-Xe 



^ (x - y- P) e~^, X e = £s + (1 - £)y, 



iPX, 



(2-4) 



where £ = corresponds to being the heavy quark coordinate while £ = 1 the light quark 
one. 

When calculating the transition amplitudes for £ = and 1, one obtains, 



J d A x (P',k\C int \P) 

= (2n) 4 5 4 (P-P'-k) [ d 3 z tr 



V/ (0)t (z; P') O^ (0) (z; P) 



(2tt) V (P-P'-k) / d 3 z tr U' (1)t (z\ P') 0^ (1) (z; P) 



-ik-z ( 2 . 5 ) 
(2-6) 



These expressions seem to be different and contradict from each other. However, when one 
rewrites Eqs. (12- 5[ I2-6|) in terms of the rest frame wave functions, these become the same 
expression given by 



(G-'OG)^ I dVV^OT'jM) ^(z'-M) e ~ 2m ^ 



(2-7) 
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where G is a Lorentz transformation matrix. 

The rest of our formulation is based on our former paper® in which it is shown how 
to construct the Lorentz-boost wave function from the rest frame wave function. In that 
paper, it is also shown that the Breit frame only with t' = x ' = y ' , i.e., times of the final 
heavy quark Q(y ' ,y ') and light anti-quark q(x ' ,x ') are equal to each other, gives the 
Lorentz-invariant results, and hence in this paper we also adopt this frame. We study main 
decay modes of the heavy-light mesons, that is, one heavy-light meson decays into another 
with one chiral particle or one photon (radiative decay), e.g., D s (0 + ) — > D s (0~) + tt° or 
D s (0 + ) — > D s (l~) + 7. The amplitudes for each process can be written in terms of form 
factors which are given in the next section. 



§3. Form Factors 



In this section, we study the form factors with an initial state of spin up to 2 + . In the 
case of radiative decays, amplitudes are given in terms of thirty electromagnetic form factors 
as follows (if one takes into account the gauge invariance, the number of independent form 
factors is reduced to twenty-three.) : 
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(3-1) 
(3-2) 
(3-3) 

(3-4) 



(3-5) 
(3-6) 



(3-7) 
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where = (— e q )q c 'y ll q c + cqQ^^Q with e q and tQ being electric charges of q and Q, re- 
spectively. In the case of decays with one chiral particle, amplitudes are given in terms of 
twenty-six chiral particle form factors as follows: 



(Q~bV|i-) 

(o-b5 f jo + ) 

i y /M 2 M 1 

<i-|j5m|Q + ) 
<Q-|j5m|i + > 

y/M 2 Mi 
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(3-8) 

(3.9) 

(3-10) 
(3-11) 

(3-12) 
(3-13) 



.(9) 
=A2 



1 aA */9 7 5 

+e a/ 3 7 , 5 e 1 V2\e 2 H v{v 2 



(3-14) 



Here and e M are polarization vectors, which could have subindex to distinguish whether it 
is for initial [i = 1) or final (2 = 2). All form factors depend only onw = (v 1 -t^), where i>j are 
initial and final velocities of heavy-light mesons, Hi (Hi — > H 2 + tt /7). These amplitudes can 
be written in terms of form factors a la semi-leptonic decay of a heavy-light meson. These 
are, however, not described in terms of just one form factor like the Isgur-Wise function 
because the interaction can not be regarded as a point contrary to the semi-leptonic decay. 

All amplitudes given above Eqs. fl3-ll) - f)3-14p are in principle calculable at the lowest order 
in 1 / rriQ because we know the explicit forms of wave functions of any spin in the heavy quark 
symmetric 

limit™™ w hi c h are described in the next section. 

54. Sum Rule 



In this section, we show one of our results for the radiative decay of a heavy-light meson. 
By analytically calculating amplitudes given by Eqs. fl3Tp - fl3-7l) . we have the following sum 
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rule in the heavy quark symmetric limit: 

p ( " 3 P!" (i+) -> l s (o-) + 7) + r ( " 3 p x " (i + ) ->• 3 Si (1-) + 7 ) 



^(0+)^ (i-)+7) = ^ 



r ( -ipx" (i + ) % (0-) + 7) + r ( "^i" (i+) 3 5 x (1-) + 7 ) 



p( 3 p 2 (2+)^ 3 5i (i-)+ 7 ) = h 



2e 2 fc 3 



9 2tt 



y rfr (« 2 1 tt 1 2 + w 2 \ 2 ) • 



(4-1) 



(4-2) 



In the heavy quark symmetric Hamiltonian, the light quark degrees of freedom are conserved 
and the LS q couling is included. Hence, the mass eigenstates in the heavy quark symmetric 
limit become mixed states of pure states,™! and are given by 






(4-3) 



s/2 -1 
1 V2 

In Eqs. f l4-l[l4-2p . u k (r) and vf(r) are initial (i = 1) and final (i = 2) radial upper/lower parts 
of the wave functions with a quantum number k as explained below. The wave function, 
which has two spinor indices and is expressed as 4 x 2 matrix form, is explicitly given by 

u k (r) 



V k (r) 

7 m\ I 



2M1 



aj mi 



(4-4) 



An r \ -iv k (r) (a ■ n) 

where n — r/r, y k m being of 2 x 2 form are the angular part of the wave function, and u k (r) 
and v k (r) are the radial parts. The total angular momentum of a heavy meson J is the sum 
of the total angular momentum of the light quark degrees of freedom j q and the heavy quark 
spin \E Q : 

1 

lq ' vv±u±± Jq ~ L ^ 2 

where \S t ' 1 

light antiquark, respectively. Furthermore, k is the quantum number of the following spinor 
operator j^HS 



L + -E q , 



(4-5) 



• -"/ ■ " 2^9 -^2x2) and L are the 4x4 spin and the orbital angular momentum of a 



K 



k\P k 



(4-6) 



j q and k are good quantum numbers in the heavy quark limit, hence the wave function can 
be described in these two. 

Applying Eqs. (14- ip and (14 -2p to actual heavy-light mesons, e.g., to D s , we have 

P (P sl (2460) ->• P s (1968) + 7) + P (P sl (2460) ->• P; i (2112) + 7) 

= P (P s0 (2317) -> D; 1 (2112) + 7) , (4-7) 

P (P si (2536) -> P s (1968) + 7) + P (P si (2536) -> £> s * 1 (2112) + 7) 

= P (P s2 (2573) -> D; 1 (2112) + 7) , (4-8) 
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where experiments give 

r (£> sl (2460) ->■ £> s (1968) + 7) < 630 keV, (4-9) 

r (£> sl (2460) -»■ ^(2112) + 7) < 280 keV, (4-10) 

r (L> s0 (2317) D*!(2112) + 7) not yet seen, (4-11) 

r (£> sl (2536) -»■ D s (1968) + 7) not yet seen, (4-12) 

r (£> sl (2536) ->■ D* 1 (2112) + 7) possibly seen, (4-13) 

r (£> s0 (2573) -> D* 1 (2112) + 7) not yet seen, (4-14) 



We have to wait and see to what extent the sum rules Eqs. (14- 1[ 14-21) are satisfied until 
experiments observe the above and other radiative decay modes. 

§5. Summary 

We have proposed a new method given by Eq. (I2-3P how to calculate decay rates for a 
heavy-light mesons and have applied it to the processes H\ — » H2 + <p a /l with one chiral 
particle <p a or one photon in the final state. We have derived form factors of vector as well 
as axial vector currents, which can be used to calculate the above processes. 

We have derived the sum rules given by Eqs. (14-11) and (14-21) for the radiative decay rates 
in the heavy quark symmetric limit, 

J2 r( u3 p^ (i + ) ^5 + 7) = r(o+ ->■ 1- + 7) , (5-i) 

5=0-, 1- 

r ( ia Pi" (i+) s + 7) = r (2 + -> r + 7 ) , (5-2) 

s=o-,i- 

which needs to be checked by future experiments. We still need to give numerical results for 
all the processes in concern, H\ — > H2 + a /7, which will be published in near future.^ 
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